X-ray flares following short 7-ray bursts from shock 
heating of binary stellar companions 

Andrew I. MacFadyen^, Enrico Ramirez-Ruiz^''^ & Weiqun Zhang^'^ 

^Institute for Advanced Study, School of Natural Sciences, Einstein Drive, Princeton, N.J. 08540 
^Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, P.O. Box 20450, 
MS 29, Stanford, CA 94309 
^Chandra Fellow 

The discovery of long-lasting (~ 100 s) X-ray flares following short gamma-ray bursts ini- 
tially called into question whether they were truly classical short-hard bursts^'^. Opinion 
over the last few years has coalesced around the view that the short-hard bursts arise from 
the merger of pairs of neutron stars, or a neutron star merging with a stellar-mass black 
hole^"^. The natural timescales associated with these processes*, however, essentially pre- 
clude an X-ray flare lasting ~ 100 s. Here we show that an interaction between the GRB 
outflow and a non-compact stellar companion at a distance of ~ a light-minute provides a 
natural explanation for the flares. In the model, the burst is triggered by the collapse of a 
neutron star after accreting matter from the companion. This is reminiscent of type la su- 
pernovae, where there is a wide distribution of delay times between formation and explosion, 
leading to an association with both star-forming galaxies and old ellipticals. 

Until recently, short GRBs were known predominantly as bursts of 7-rays, largely devoid of 
any observable traces at any other wavelengths. However, a striking development in the last several 
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months has been the measurement and locaUzation of fading X-ray signals from several short 
GRBs, making possible the optical and radio detection of afterglows, which in turn enabled the 
identification of host galaxies at cosmological distances^~^°. The presence in old stellar populations 
e.g., of an elliptical galaxy for GRB050724, rules out a source uniquely associated with recent 
star formation^ ^ In addition, no bright supernova is observed to accompany short GRBs^'^'^^, 
in distinction from most nearby long-duration GRBs^^. The current view of most researchers is 
that short GRBs arise when a neutron star (NS) binary or neutron star-black hole (BH) binary, 
which loses orbital angular momentum by gravitational wave emission, undergoes a merger^"^. 
Current calculations of compact binary mergers suggest that high spatial velocities would take 
these binaries, in more than half of the cases, outside of the confines of the host galaxy before they 
merge and produce a burst, in agreement with current observations^"*' 

Recently, long duration (~ 100 s) X-ray flares have been observed to follow several short 
GRBs*^'*^, e.g., GRB050724, after a delay of ~ 30 s. There is also independent support that 
X-ray emission on these timescales is detected when lightcurves of many bursts are stacked*^. 
One possible interpretation of these (rapidly declining) flares is that a large fraction of energy 
continues to be emitted by the GRB source for as long as minutes*^. This hints at the desirability 
for a "central engine" lasting much longer than a typical dynamical time scale for a stellar mass 
compact object. This is in disagreement with theoretical estimates, which suggest that NS-NS 
and NS-BH mergers will lead to brief energy input episodes, typically of the order of the duration 
of a short burst^. It is argued here that the flare emission is, however, naturally produced by the 
interaction of the extended - possibly magnetically dominated - ejecta emanating from a short 
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GRB with the envelope of a giant companion star. This more isotropic GRB ejecta component 
need not necessarily dominate the total burst energetics, but it can be efficiently reprocessed by 
the envelope of the companion star at distances ~ 10^^ cm, into an X-ray flare with a luminosity 
and timescale comparable to the observed values. The portion of the ejecta shell not impacting 
the star expands undisturbed to large radius (10^'' — 10^^ cm) where internal shocks or magnetic 
dissipation convert its energy into 7-ray photons as envisioned in standard models^^. The observer 
is envisioned here to be located out of the plane of the binary. Soft X-rays produced during the 
collision, although produced at smaller radii than the 7-rays, are observed to arrive after the GRB. 
Under this interpretation, the GRB engine is not required to operate for longer than a typical 
burst duration of 0.1-1 s. We note here that, in addition to the rapidly declining ~ 100s flare, 
GRB 050724 rebrightens^^ after ~ 10^ s. Given its slower decline, this feature has a satisfactory 
standard^^ explanation: as the leading edge of fast ejecta moves farther away from the central 
engine, it starts to decelerate and is caught up by slower-moving material, creating a 'refreshed 
shock' . 

Three snapshots from numerical simulations of a GRB ejecta shell as it expands, makes con- 
tact with a red giant star, and ultimately engulfs it, are shown in Figure 1. The interaction begins 
with the formation of a strong shock as the shell is rapidly decelerated from Fj ~ 400 to F/ ~ 1 
in the surface layers of the star, dissipating a large fraction of its kinetic energy. Internal energy is 
created rapidly and remains at a few 10^^ erg for at least ~ 100 s (Figure 2). The initial rise in in- 
ternal energy takes place immediately after the shell reaches the stellar surface, 30 s after the GRB. 
A strong shock forms, starting where the two spheres first touch, then subsequently spreads as an 
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increasingly larger portion of the stellar surface is hit. The rate of energy dissipation as the shell 
sweeps across the star - for a GRB located approximately one stellar radius from the stellar surface 
- is roughly given by E Ec/{AR^) = 7.5 x 10^^(E/10^^ erg)(i?*/10^2 cm)-^ erg/s, where 
E is the total energy in the ejected shell, c is the speed of light and is the stellar radius. This 
simple estimate yields luminosities that are consistent with those found in numerical calculations 
shown in Figure 2. A few hundred seconds after the GRB (bottom panel of Figure 1), the shell 
has almost entirely wrapped around the star and fresh material is no longer being shocked. The 
internal energy declines rapidly after the GRB flare finishes shocking the star bringing the flare 
emission to an abrupt end. The decay slope of the corresponding lightcurve may be substantially 
modified by radiative effects. However, the dissipated energy has the correct magnitude, delay and 
duration to account for the observed properties of X-ray flares in short GRBs. 

The flares observed after the GRB 050724 and GRB 050709 bursts corresponds, for their 
assumed distances of z — 0.257 and z — 0.16, to £■/ ~ 6 x 10^^ erg and Ef r-^ 3 x 10^^ erg 
respectively in the 2 — 25 keV band. How is energy dissipated during the GRB-star interaction 
(Figure 2) transformed into soft X-ray radiation? One possibility is that the GRB ejecta shell de- 
velops a stand-off shock before encountering the stellar envelope. The post-shock region generates 
turbulent magnetic fields and accelerates electrons which produce a synchrotron power-law radi- 
ation spectrum. The magnetic field strength would be of order 10^ G at 10^^ cm, strong enough 
to ensure that the shock-accelerated electrons cool promptly, yielding a power-law continuum ex- 
tending into the X-ray band. Some of these X-rays would be deflected along the stellar surface 
before escaping, but about half (the exact proportion depending on the geometry and flow pattern) 
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would irradiate the material in the stellar envelope. For the high radiative efficiencies expected 
in relativistic shocks, radiative heating of the shocked material would be comparable to that of 
bulk heating^*'. However, radiative heating deposits energy near the stellar surface in layers with 
modest scattering optical depth, the temperature being determined by photoionization equilibrium. 
This shallow radiatively heated layer would be expected to be substantially hotter than a deeper 
bulk-heated region. For low radiative efficiencies, on the other hand, energy deposition by bulk 
heating would spread over a highly optically thick layer. In this case, the cooling rate, mainly 
due to bremsstrahlung, recombination and Comptonization, would be high enough to reduce the 
temperature of the bulk-heated electrons to the equivalent black body temperatures of T ~ a few 
tens of keV. This suggests that, for the conditions envisaged here, most of the flare luminosity in 
the X-ray band could be thermal. This is consistent with numerical results where a pressure of 
10^^ dyne cm~^ behind the shock (see the middle panel of Figure 1), corresponds to a black-body 
temperature of ~ 2.5 x 10^ K and a characteristic photon energy of ~a few keV. We note that a 
radiatively-heated layer with a density up to ~ 10^^ cm~^, could produce strong Fe line emis- 
sion provided that the ionization parameter { = (3L/{r'^ne) exceeds^^ 10^, where L is the total 
luminosity of the GRB outflow and P is the fraction of the power that goes into X-ray continuum. 
This condition is indeed satisfied unless P < 10~^. Thermal X-ray emission could also display 
line features, and such signatures should certainly be looked for. We note that the GRB shell must 
be optically thin when its radius is ~ 10^^ cm so that the shocked star is visible through the GRB 
shell. This requires that the GRB ejecta contain 10~^ Mq or less and suggests that a magnetically 
dominated flow may be preferred^^. 
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Shock heating of a binary stellar companion has other interesting consequences. Before the 
passage of the GRB ejecta, ~ 10~^ Mq of the envelope of the companion star is compressed by 
the shock and heated. The deposited energy, a few 10^^ erg, will cause the outer layers of the star 
to expand explosively. This will not, however, produce a supernova, for two reasons. First, the 
shock temperatures are too low for radioactive elements such as Ni^^ to be produced. Second, the 
amount of ejected debris is small. Instead, after the expanded envelope becomes optically thin, a 
faint infrared/optical transient would appear a few weeks after the GRB. 

Finally, we offer a few comments on the progenitor systems in which this interaction can 
naturally occur. Models in which the short GRB results from the collapse of a rapidly rotating 
neutron star in a close binary system provide a natural scenario. In the model, the neutron star 
accretes matter from the stellar companion, eventually collapsing to form a black hole. The angular 
momentum in the equatorial region of the rapidly-spinning system is too large to be swallowed 
immediately when the black hole forms. The expected outcome, after a few milliseconds, would 
therefore be a spinning black hole orbited by a torus of neutron-density matter. The mass in these 
disks is found to range from^-' ~ 10~^ — 1Q~'^Mq. If magnetic fields anchored in the disk 
do not thread the black hole, then a relativistic outflow powered by torus accretion can at most 
carry the gravitational binding energy of the torus. The extractable energy in this case is several 
times 10^° e{Mi/ 10^^ Mq) ergs, where e is the efficiency in converting gravitational energy into 
relativistic outflow. If magnetic fields of comparable strength thread the black hole, its rotational 
energy offers an extra (and even larger) source of energy that can in principle be extracted via the 
Blandford-Znajek mechanism^'*. 
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Not surprisingly, there is more than one way to produce a rapidly rotating, neutron star in 
a close binary including, for example, a common envelope evolution in which a neutron star 
is enveloped by the expanding atmosphere of a giant companion. Alternatively, a binary initially 
containing two massive stars, could form a neutron star and helium star binary. During the common 
envelope phase, the neutron star may accrete over one solar mass and coUapse^^. Thus, we expect 
these systems to be associated with star forming regions. Reminiscent of what is observed in 
type la supemovae, the model predicts a wide distribution of delay times between formation and 
explosion, and in turn the association of short hard GRBs with both star forming galaxies and with 
ellipticals dominated by old stellar populations^''. 

Neutron stars inspiralling into a stellar envelope can accrete at rates vastly exceeding the 
Eddington limit if the flow develops temperatures high enough to allow neutrinos to radiate the 
gravitational binding energy. The fate of the neutron star depends on a sensitive balance between 
the rate at which it accretes and the rate at which energy is deposited into the common envelope^^. 
Its chances for survival are therefore diminished both by increasing the accretion rate and by aug- 
menting the epoch of common-envelope evolution. Observationally, several binary pulsars are 
known whose properties are consistent with the neutron star having survived a phase of common- 
envelope evolution. Camilo et al.^' identify four pulsars which have relatively large companion 
masses in excess of O.45M0. These systems are likely to have undergone deep common-envelope 
evolution with low-mass companions'^ (1 — 3Mq). If the neutron star is to be able to accrete a 
large mass during inspiral, the initial distance between the two stars after the first mass transfer 
exchange should not be much larger'^ than R^. This ensures that a companion star struck by the 
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expanding GRB outflow will intercept a large fraction of the ejected material, thus ensuring the 
production of a bright flare. 

Much progress has been made in understanding the nature of cosmological gamma-ray 
bursts. Still, various alternative ways of triggering the explosions responsible for short GRBs 
remain: NS-NS or NS-BH binary mergers, spun-down supramassive NS^*^ and accretion induced 
collapse of a NS. The presence of X-ray flares may help distinguish between viable progenitors. 
In the absence of a supernova-like feature, the interaction of GRB ejecta with a stellar binary com- 
panion may be the only observable signature in the foreseeable future shedding light on the identity 
of the progenitors of cosmological short bursts. 
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Figure 1 



Collision of GRB ejecta with a companion red giant star showing the logarithm of pressure 
at t = 10s (top), t = 42s (middle) and t = 184s after the GRB . The observer is located envisionedo 
be out of the binary plane toward the top of the figure. In the top panel, the small orange circle 
is the GRB ejecta and the larger purple sphere is a red giant star with radius 1.5 x 10^^ cm. The 
ejecta form a spherical shell containing 2 x 10^^ g (10~^Mq) expanding at nearly the speed of 
light. 30 seconds after the GRB, the ejecta collide with the surface of the red giant and are rapidly 
decelerated in a strong shock. The internal energy produced during the collision is capable of 
explaining the energetics of the observed X-ray flares. The middle panel shows the GRB blastwave 
and shocked star 12 seconds after impact. The dark red region is the X-ray emitting region. 
In the bottom panel, the blastwave has wrapped around the star. At this time, the energy of 
the blastwave is mainly kinetic. The orange region in the middle panel and the green region in 
the bottom panel show the reflected ejecta shell traveling leftward toward to GRB source. The 
undisturbed shell expands close to the speed of light for at least two orders of magnitude in radius 
before producing the prompt GRB. We used the RAM code^° to solve the conservative equations 
of special relativistic hydrodynamics. The computational domain for the simulations was < 
r < 1.2 X 10^^ cm with a maximal resolution of Ar = 4.7 x 10''' cm for 1 D and -1.2 x 10^^ 
cm< z < 1.2 X 10^3 cm, < r < 1.2 x 10^3 cm, with Az = Ar = 7.32 x 10^ cm for 2D. A 
gamma law equation of state with r = 4/3 was used and the ambient density was 100 baryons 
cm~3. 

Figure 2 Total internal energy. As the cold GRB shell collides with the stellar surface a large 
fraction of the energy in the blast wave is converted to internal energy in a shock. Internal energy 
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is produced with magnitude, duration and delay appropriate for being tine source of X-ray flares 
observed to follow short GRBs. The sharp decline in the flare is easily explained by the finite time 
the companion's surface is actively shocked. In this figure, the remaining internal energy is being 
reconverted to kinetic energy and is not expected to contribute to the flare lightcurve. 
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